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Abstract. The distance to NGC 5128, the central galaxy of the Centaurus group and the nearest giant elliptical 
to us, has been determined using two independent distance indicators: the Mira period-luminosity (PL) relation 
and the luminosity of the tip of the red giant branch (RGB). The data were taken at two different locations in 
the halo of NGC 5128 with the ISAAC near-IR array on ESO VLT. From more than 20 hours of observations 
with ISAAC a very deep its-band luminosity function was constructed. The tip of the RGB is detected at 
K s = 21.24 ± 0.05 mag. Using an empirical calibration of the it-band RGB tip magnitude, and assuming a 
mean metallicity of [M/H] — —0.4 dex and reddening of E(B — V) = 0.11, a distance modulus of NGC 5128 of 
(m — M)o = 27.87 ± 0.16 was derived. The comparison of the ii-band RGB tip magnitude in NGC 5128 and 
the Galactic Bulge implies a distance modulus of NGC 5128 of (m — M)o = 27.9 ± 0.2 in good agreement with 
the it-band RGB tip measurement. The inner halo field has larger photometric errors, brighter completeness 
limits and a larger number of blends. Thus the RGB tip feature is not as sharp as in the outer halo field. The 
population of stars above the tip of the RGB amounts to 2176 stars in the outer halo field (Field 1) and 6072 
stars in the inner halo field (Field 2). The large majority of these sources belong to the asymptotic giant branch 
(AGB) population in NGC 5128 with numerous long period variables. Mira variables were used to determine 
the distance of NGC 5128 from a period-luminosity relation calibrated using the Hipparcos parallaxes and LMC 
Mira period-luminosity relation in the it-band. This is the first Mira period-luminosity relation outside the Local 
Group. A distance modulus of 27.96 ±0.11 was derived, adopting the LMC distance modulus of 18.50 ±0.04. The 
mean of the two methods yields a distance modulus to NGC 5128 of 27.92±0.19 corresponding to D — 3.84±0.35 
Mpc. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: stellar content - Stars: fundamental parameters - 
Galaxies: individual: NGC 5128 



| 1. Introduction 

Thanks to its proximity NGC 5128 (Centaurus A) has 
attracted lots of attention. It is the closest representative 
of the active radio galaxy and giant elliptical galaxy class 
of objects and is among the nearest AGNs. The latest 
review (Israel IT~998p summarises its characteristics and it 
is a good starting point for the vast literature about this 
galaxy. It also gives a summary of distance determinations, 
which started converging in the last decade from the early 
highly uncertain values ranging between 2.1 and 8.5 Mpc 
(Sersic IT958I Sandage & Tammann I1974f) towards a less 
uncertain range between 3.2 Mpc (Hui et al. I1993f) and 
4.2 Mpc (Tonry et al. l2UtIT)l . 
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* Based on observations collected at the European 
Southern Observatory, Paranal, Chile, within the Observing 
Programmes 63.N-0229, 65.N-0164, 67.N-503, 68.B-0129 and 
69.B-0292 and at La Silla Observatory, Chile, within the 
Observing Programme 64.N-0176(B). 



The recent determinations of distance to NGC 5128 
applied a range of methods. Hui et al. I|1993l used the 
planetary nebula luminosity function to derive a distance 
modulus of (m - M) = 27.73 ± 0.14 mag. The globu- 
lar cluster luminosity function was analysed by Harris et 
al. 1|1988[> yielding (m-M) = 27.53±0.5 mag. More recent 
globular cluster searches in this galaxy (e.g. E.eikuba l2()()ll 
Peng 2003 ) will allow a more precise determination of dis- 
tance with this method through a much better sampled 
globular cluster luminosity function. 

The luminosity of the red giant branch (RGB) tip stars 
in the /-band is a recognised distance indicator (e.g. Lee 
et al.Ei23- In NGC 5128 it was first used by Soria et al. 
ljTTJ9l)|) who resolved the stellar halo using HST±WFPC2. 
They derived a distance modulus of (m — M) = 27.86 ± 
0.16 mag for WF chips and (m - M) = 27.76 ± 0.16 mag 
for the PC chip, and adopted a mean distance modulus 
of (m — M)o = 27.8 ± 0.2 mag. More recently, Harris et 
al. fiMty used deeper HST±WFPC2 photometry in a less 
crowded field to derive the distance. Their RGB tip lumi- 
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nosity analysis resulted in (m — M)o = 27.98 ± 0.15 mag. 
Moreover, the same authors adjusted the Hui et al. (1993) 
distance modulus to (m — M) = 27.97 ± 0.14, increas- 
ing it by 0.2 mag in order to correct for the contempo- 
rary Local Group distance scale and the M31 distance of 
(m - M)m31 = 24.5 (van den Bergh lTMSl Harris H^|l . 

More than the planetary nebula luminosity function, 
the surface brightness fluctuations (SBF) method received 
several revisions of its calibration zero point. Tonry & 
Schechter ifTWij) first derived (m - M) = 27.48 ± 0.06 
using I-band SBF. This value was subsequently revised to 
(m - M) = 27.71 ± 0.10 by Israel l(TM5|l who used the re- 
sults from Tonry lfTMT)l and then to (m - M) = 28.18 ± 
0.07 by Marleau et al. pUUUjl after Tonry et al. (|TM7|) . 
Most recently, Tonry et al. 1)2001[) report /-band SBF 
yielding a distance modulus of (m — M) = 28.12±0.15. In 
their list of nearby galaxies with SBF distance measure- 
ments, NGC 5128 occupies 9 th place and is the nearest 
giant elliptical galaxy. The distance modulus determina- 
tions from the literature are summarized in Tabled 

In this paper I use Mira variables from the long 
period variable star catalogue in NGC 5128 (Rejkuba 
et al. I2003a|) and the H and K-haxid luminosity func- 
tions to derive independent measurements of distance to 
NGC 5128. Data are briefly described in Sect. and in 
Sect. J, H, and X-band luminosity functions are anal- 
ysed. The distance to NGC 5128 is derived from the RGB 
tip magnitude in Sect. 14.11 The Mira period- luminosity 
(PL) diagram is constructed in Sect. 14.21 and used to de- 
termine the distance to NGC 5128. Finally the results are 
summarized in Sect. 

2. The data and photometry 

Data used here were described in detail in Rejkuba et al. 
(EHH1 and EIIS3I ■ They consist of a set of 20 if s -band 
epochs, 1 J s , and 1 ff-band image of a field located ~ 17' 
north-east from the center of NGC 5128 (Field 1). The 
second field (Field 2) was observed once in the J s and H- 
bands and 24 times in the -ftT s -band. It is located ~ 9' south 
from the galactic center. All the observations, except one 
i^-band epoch of Field 2, were taken with the ISAAC 
near-IR array at UT1 (Antu) telescope at ESO Paranal 
Observatory in service mode. The additional if s -band ob- 
servation of Field 2 was secured during an observing run 
with the SOFI near-IR array on the NTT at ESO La Silla 
Observatory under exceptional seeing conditions. 

For the details about the photometry and completeness 
simulations the reader is referred to the above mentioned 
papers. The photometry is complete more than 50% in the 
K s and iJ-bands for stars brighter than 22.5 for Field 1 
and 21.6 for Field 2. The 50% completeness limit for the 
J s -band is 23.25 and 22.5 mag for Field 1 and 2, respec- 
tively. 

Multi-epoch if s -band observations span a 1197 day in- 
terval and were used to search for variable stars. A total 
of 1504 red variables were detected in the two halo fields. 
For 1146 variables with at least 10 good measurements, 
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Fig. 1. Color-color diagram for all the stars in Field 1 and 
2 with ALLFRAME photometry errors smaller than 0.1 
mag. Foreground Galactic dwarf and giant stars follow the 
long-dashed and solid lines indicating their intrinsic colors 
(Bessell & Brett 1988). A reddening vector corresponding 
to Fi(B — V) = 0.5 mag is shown. 

periods, amplitudes and mean K s magnitudes were de- 
termined using Fourier analysis and non-linear sine-curve 
fitting algorithms. Almost all these variables belong to 
the class of long period variables (LPVs) with Mira and 
semiregular variable stars. The complete catalogue of all 
LPVs is presented in Rejkuba et al. (2003a) to which 
reader is referred for discussion of completeness and ac- 
curacy of period and amplitude determination. 

The ISAAC filters used are J s , H and K s . The J s filter 
was preferred over J, due to the red leakage of the latter. 
It is centered at 1.24 /jm and has a width of 0.16 /xm. 
The transformation to the J-band of the LCO photomet- 
ric system (Persson et al. I1998fl was obtained through the 
following transformation (Chris Lidman private commu- 
nication): 

Jlco = Jsisaac + 0.033 * (Js - Ks)jsaac - 0.022 (1) 

The ISAAC H and K s filters are on the LCO photometric 
system. 

It should be noted that to all the if s -band magnitudes 
from the catalogue (Rejkuba et al. 2003a) a constant of 
0.1 mag has been subtracted. A comparison of the color- 
color diagrams with the intrinsic colors of Galactic dwarfs 
and giants (Bessell & Brett I1988|) has revealed an error 
in the aperture correction of the if s -band photometry. 
For reference the J — H vs. H — K color-color diagram is 
shown in Fig. with the Bessell & Brett fiducials 
overplotted. The solid line indicates the intrinsic colors of 
early to late type giants and the long-dashed line is for 
dwarfs. Foreground Galactic stars are found along these 
lines. The large majority of the stars in NGC 5128 are 
found close to the locus of late type giants and redwards 
from there, where long period variables and carbon stars 



M. Rejkuba: Distance to NGC 5128 
Table 1. Summary of the distance distance modulus (DM) determinations from literature. 



3 



# 


DM (mas) 


Method 


Reference 


1 

1 


zu.u 


>"i"^ lor 111 l~V~t l»"l*""\C?l"d""\7" TllT~lf'*"d~l/""\Tl f 1 H 1 

OLclldl lUXXl-UlUolUy ILIIILLIUII 1±j_T 1 


ucl alL lyoo 


9 


9Q 


Largest HII regions 


oaiiQatie ol iammanniiy/''±i 


o 
o 


97 70 J.A 1/1 


Planetary nebula, LF 


riui ex ai. iyyo 






Planetary nebula LF 


(3i) revised bv Harris et al. 119991 


r. 



Z ( .Do zt U.O 


Globular cluster LF 


riarris ei ai. iiyooi 


6 


27 86 zb 16 


/-band RGB tip (WF chips of WFPC2) 


Soria et al. I19S6I 


7 


27.76 ±0.16 


/-band RGB tip (PC chip of WFPC2) 


Soria et al. 115351 


8 


27.98 ±0.15 


/-band RGB tip (WFPC2) 


Harris et al. TTMft 


9 


27.48 ± 0.06 


/-band SBF 


Tonry h Schechter [M5U| 


10 


27.71 ±0.10 


/-band SBF 


(9) revised by Israel 


11 


28.18 ±0.07 


/-band SBF 


(9) revised by Marleau et al.[2000 


12 


28.12 ±0.15 


/-band SBF 


Tonry et al-ISTO!! 





Fig. 2. if s -band luminosity function for all stars (solid 
black line) in Field 1 (left) and Field 2 (right). Top: the 
dotted line histogram is for stars bluer than — (V— K s )+25 
(see Rejkuba et al. 120011 . Mostly they are foreground 
Galactic stars as well as young blue and red supergiants in 
Field 1 . The variable star luminosity function (red shaded 
histogram) is also compared with the complete luminosity 
function. Bottom: A sum of the foreground and variable 
star luminosity function is compared to the total luminos- 
ity function of the two fields. 

are located. NGC 5128 photometry has been de-reddened 
adopting E(B - V) = 0.11 (Schlegel et al. 11333)1 and the 
Cardelli et al. 1)19891) reddening law and it has been trans- 
formed to the Bessell & Brett (|1988|l JHK photometric 
system using the transformation equations available from 
the 2MASS web page 1 . 



3. Luminosity functions 

The K s -ha,nd luminosity functions for all the stars in 
Field 1 (left) and 2 (right) are shown in Fig.[2as solid line 
histograms. These are not corrected for incompleteness 



nor for Galactic foreground contamination. The number 
of foreground Milky Way stars can be obtained from the 
color-color diagram (Fig.^) or from optical-near-IR color- 
magnitude diagrams published by Rejkuba et al. I|2001[l . 



There are 194 stars with near-IR colors consistent with 
dwarf or early type giants in the Field 1 color-color dia- 
gram. Similarly, 192 sources have blue colors in the VK S 
CMD (see Rejkuba et al. l2UUl|l . They are plotted as a dot- 
ted histogram in Fig. [3 These stars are mostly foreground 
Galactic stars as well as young blue and red supergiants 
in Field 1. The variable star luminosity function is plot- 
ted as a shaded histogram. It is clear that majority of 
the variables have luminosities brighter than the RGB tip 
(see below) and are in the thermally pulsing asymptotic 
giant branch (TP-AGB) phase. In the lower panels the to- 
tal luminosity functions are compared with the luminosity 
functions of the sums of the contributions of foreground 
stars and LPVs. 

There are 2176 stars in Field 1 above the tip of the 
RGB (K s < 21.24, see next section) detected in at least 
3 .ft's-band images. 198 of these sources have not been 
detected in J s and/or _ff-band frames. According to the 
VK S CMD, 192 stars have blue colors and they belong ei- 
ther to the young population in Field 1 or to the Galactic 
foreground. Among the remaining 1984 stars 426 are vari- 
able stars (accounting for 71% of variables detected in this 
field). 

In Field 2 there are 6072 sources brighter than K s = 
21.24 detected in at least 3 -RT-band images. 351 of these 
sources have not been detected in J and/or iJ-band 
frames. Of these 217 objects have blue colors and belong to 
the Galactic foreground contamination. Of the remaining 
5855 stars brighter than the RGB tip 829 are asymptotic 
giant branch (AGB) variables (accounting for 92% of de- 
tected variables in this field). 

Some of the bright non-variable AGB giants could be 
blends of 2 RGB tip stars. Renzini fI998J provided the 
equations that enable the calculation of the expected num- 
ber of stars in a given evolutionary phase (Nj) in an image 
in which the total sampled luminosity is Lt- 



1 http : //www. astro . caltech. edu/~jmc/2mass/v3/transformAtjioii^ t)Lxtj 
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Fig. 3. J s , H and K s -b&nd luminosity functions for Field 1 are shown in the top panels. Dotted lines are used for 
the observed and solid lines for the completeness-corrected luminosity functions. In the bottom panels the Sobel edge- 
detection filter response curves are displayed. Vertical dotted lines indicate 50% completeness limits and the short 
thick marks indicate RGB tip magnitudes measured in this Field. 



B(t) is the specific evolutionary flux of the population or 
the number of stars entering or leaving any post main 
sequence evolutionary stage per year and per solar lumi- 
nosity of the population. It is a very weak function of 
age and ranges from ~ 0.5 x 10~ n to ~ 2.2 x 10~ n 
stars Lq 1 yr _1 , as age increases from 10 Myr to 15 Gyr. 
tj is the duration of the evolutionary phase. For a solar- 
metallicity, 15 Gyr old population Lt — 0.36Lk- We 
can estimate the total sampled luminosity in the two 
fields by summing the observed counts on the dark sub- 
tracted, flat-fielded, sky-subtracted and calibrated frames. 
Using the combined frames for the reference epoch ob- 
served on JD=2451734, the total sampled luminosity is 
5 x 10 7 L© and 7 x 10 7 L© for Field 1 and 2, respectively. 
From the surface brightness profile of Dufour et al. I|1979J) 
the total luminosities are 5.9 and 28.7 x 10 7 L©, while fr om 
that of Mathieu et al. (fTMfJ|) they are 9.5 and 31 x 10 7 LQ. 
Our measured value for Field 1 is close to that of Dufour et 
al., but the measured luminosity of Field 2 is much lower. 

From Table 1 of Renzini ( 1998} the expected num- 
ber of LPVs in an old near-solar metallicity population 
is ~ 5/10 6 LQ. The total sampled luminosity of Field 1 
is 50-100 times higher, yielding 250-500 LPVs in agree- 



ment with the 437 LPVs detected in that field (Rejkuba 
et al. |2DP3a) . In Field 2 our measured luminosity would 
suggest ~ 350 LPVs, but the surface brightness measure- 
ments from the literature predict ~ 1500 LPVs. For com- 
parison, we have detected 903 variable stars and measured 
periods for 709 LPVs in Field 2. 

The probability that a pixel contains two stars in phase 
7 is ~ TV? and the number of 2-star blends in a frame with 
a total number of pixels N plx is: 

N 2j = N*N pix = [B(t)LTtj} 2 /N pix (3) 

The resolution of ground based images is worse than 1 
pixel and the sampled luminosity has to be divided by the 
available number of resolution elements. Each resolution 
element is 0'.'31 for Field 1 and 0'.'36 for Field 2, corre- 
sponding to the best seeing epochs used to detect sources. 
Hence the expected number of blends of two RGB tip 
giants in Field 1 is between 100 and 400 for the total sam- 
pled luminosity between 5 and 9.5 x 10 7 LQ. In Field 2 it 
is ~ 250 if the measured luminosity from the image is as- 
sumed, but it rises to 4900 for 31 x 10 7 LQ from Mathieu 
et al. The expected number of blended stars in the vicin- 
ity of the RGB tip is low enough to measure the RGB 
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Fig. 4. J s , H and X s -band luminosity functions for Field 2 are shown in the top panels. Dotted lines are used for 
the observed and solid lines for the completeness-corrected luminosity functions. In the bottom panels the Sobel edge- 
detection filter response curves are displayed. Vertical dotted lines indicate 50% completeness limits. Short thick marks 
indicate RGB tip magnitudes measured in Field 1 (Fig. 01. 



tip discontinuity in the luminosity function of Field 1, but 
may smear out that feature in Field 2. Taking the highest 
probable number of blends, there are at least 1150 non- 
variable bright AGB giants in Field 1 and some 150 in 
Field 2. Rejkuba et al. l|2003ajl made extensive simula- 
tions and estimated the completeness of the detection of 
LPVs as a function of their magnitude, period and ampli- 
tude. These simulations show a similar completeness level 
of the LPV catalogue in the two fields implying that the 
number of bright extended giant branch non- variable stars 
is larger in Field 1. However, these stars could be semi- 
regular small-amplitude (AK ^ 0.3) variables for which 
the LPV catalogue is less complete. Their near-IR prop- 
erties will be discussed together with near-IR properties 
of the LPVs in a forthcoming paper (Rejkuba et al. in 
preparation). 



LPVs account for at least 26% and 70% of the extended 
giant branch population in Fields 1 and 2, respectively. 
Their high luminosity cut-off extends beyond Mk = —8.7. 
The maximum brightness achieved by AGB stars is evi- 
dence for the presence of an intermediate-age population. 



4. NGC 5128 distance 

In this section the distance to NGC 5128 is measured using 
two different methods: (i) the brightness of the RGB tip 
in the K and iJ-bands and (ii) the Mira period-luminosity 
relation. 

4.1. RGB tip 

The J s , H and i^-band luminosity functions of Field 1 
stars are presented in Fig. 01 These are the deepest lu- 
minosity functions for a stellar population belonging to a 
galaxy beyond the Local Group. The solid lines are used 
for completeness-corrected luminosity functions and dot- 
ted lines for the observed ones. Vertical dotted line in- 
dicates 50% completeness limit of the photometry. Short 
thick marks indicate the measured RGB tip magnitudes 
in Field 1. 

In the bottom panels of Fig. [3] I show the luminosity 
functions convolved with the Sobel edge detection filter 
with a (-2, 0, +2) kernel (Lee et al. 1199311 . They are used 
as localised slope estimators with 2-point smoothing. 
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There is a very sharp peak in the Field 1 response 
function curve at J = 22.54 ±0.06, H = 21.54 ±0.06, and 
K s = 21.24 ± 0.05. For easier comparison with the models 
and literature, the J-band measurement is transformed to 
the LCO photometric system according to Eqn. ^ These 
maxima appear at the position where there is a strong 
change in the slope of the RGB luminosity function and it 
corresponds to the magnitude of the tip of the RGB. The 
error in the RGB tip magnitude is obtained by averag- 
ing over 100 different measurements of Sobel kernel peak 
magnitude, where luminosity function sampling is varied 
by 0.001 mag over 0.08 to 0.2 mag. 

The RGB tip feature is smoothed over a couple of 
tenths of magnitude in Field 2 (Fig. 0J) . From the 2-point 
Sobel kernel, the RGB tip in Field 2 is centered around 
K s = 21.22 ± 0.15 mag. The broadening and smoothing 
of the peak is due to much shallower completeness lim- 
its (note the short thick marks indicating the RGB tip 
magnitude measured in Field 1), resulting in larger pho- 
tometric errors, as well as due to contribution of blends, 
the presence of AGB stars and a larger metallicity spread 
in this field. 

The tip of the RGB in the iif-band is a function of 
metallicity. Its dependence has been measured empirically 
by Ferraro et al. ( 2000) using a set of Galactic globular 
clusters with a range of metallicities. They obtained the 
following relation: 



M 



TRGB 



= -(0.64 ± 0.12) [M/H] - (6.93 ± 0.14) 



(4) 



The mean metallicity for the stars in the NGC 5128 
halo has been inferred to be around [M/H] = —0.4 dex 
(Marleau et al. 12170171 Harris & Harris 12170171 Rejkuba 
2002). With this estimate, and assuming a reddening of 
E(B - V) = 0.11 (Schlegel et al. IT55g)l . corresponding 
to A K = 0.039 (Cardelli et al. I1989[> . the empirical cali- 
bration of the RGB tip i^T-band magnitude can be used 
to derive the distance to NGC 5128. The difference be- 
tween the photometric bands is of the order of 0.01 mag 
(Persson et al. I1998f) . This yields a distance modulus of 
(m — M)o = 27.87 ± 0.16. The quoted error does not in- 
clude the uncertainty in reddening and in mean metallicity 
of the stars. A higher mean metallicity by 0.1 dex implies 
a larger distance modulus by 0.06 mag. The -RT-band is 
rather insensitive to extinction, and reddening as high as 
E(B — V) = 0.35 would lower the distance modulus by 
0.088 mag. 

The above measurements can be compared to that of 
the RGB tip brightness in the LMC and in the Galactic 
Bulge. For example, Cioni et al. (20001 have measured 
K s rRGB (LMC) = 11.94 ± 0.04 and J TRGB {LMC) = 
13.06 ±0.02; these values are not corrected for metallicity. 
The difference between the RGB tip magnitudes in the 
LMC and NGC 5128 imply distance moduli differences 
of A(m - M) Kg = 9.26 ± 0.06 mag and A(m - M)j = 
9.39 ± 0.06 mag. Adopting a distance modulus of 18.5 for 
the LMC (Alves et al. 12002(1 yields distance moduli for 
NGC 5128 of 27.76 and 27.89 using K s and J bands, re- 
spectively. While this is within errors consistent with the 




Fig. 5. Period-luminosity diagram for all variables in 
NGC 5128 with 10 or more data points with well de- 
termined periods from the LPV catalogue of Rejkuba et 
al. I|2003a[) . Field 1 variables are plotted with filled and 
variables in Field 2 with open symbols. 



above quoted 27.87 ± 0.16, it should be remembered that 
the higher metallicity of the NGC 5128 stars with respect 
to the LMC, makes the RGB tip magnitude brighter in K- 
band (e.g. Zoccali et al. l2003f) . hence the relative difference 
of distance modulus is a lower limit in the if s -band. 

The J-band RGB tip magnitude is expected to de- 
crease slowly as the metallicity increases, while the tip 
luminosity in the if-band appears to be fairly constant 
in the metallicity range —1.0 [Fe/H] ^ 0.0 (Zoccali 
et al. I2003J) . Hence the latter could be used for distance 
determination. Schulte-Ladbeck et al. I|1999|l have investi- 
gated the behaviour of the ii-band tip luminosity with age 
and metallicity. They recommend Mjj RGB = —5.5 ±0.1 
for -2.3 < [Fe/H] < -1.5. Bertelli et al. {HSU 



models 

suggest Mjj RGB ~ — 6 mag for higher metallicity, while 
m trgb = _ g 4 ± q 2 i n the Galactic Bulge (Zoccali 
et al. 2003). The latter yields a distance modulus to 
NGC 5128 of 27.9 ± 0.2. Due to high extinction, small 
number of stars and possible differential reddening in 
the Galactic Bulge field observed, the M]j RGB value is 
rather uncertain. However, it yields a distance modulus 
to NGC 5128 of 27.9 ± 0.2 in excellent agreement with 
that derived from the if-band RGB tip magnitude. 

4.2. Mira PL relation 

In Fig. 01 plot the mean magnitudes vs. logarithm of pe- 
riods for all the variables with 10 or more data points for 
which derived periods had a significance parameter from 
Fourier analysis < 0.7 (Rejkuba et al. 2003a). Field 1 vari- 
ables are shown with filled and Field 2 with open sym- 
bols. The mean magnitudes are corrected for extinction 
and transformed to the SAAO photometric system (Carter 
1990) in which most Mira studies are published. 
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Fig. 6. Example of two LPVs where aliasing is present. 
Both periods, a short one of ~ 1/2 yr and ~ twice as long 
one can be equally well fit. 
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Fig. 7. Period-luminosity diagram for all variables in 
NGC 5128 with 10 or more data points with well deter- 
mined periods. Alternative aliasing periods as listed in 
Table |3 are plotted. Field 1 variables are plotted with 
filled and variables in Field 2 with open symbols. 



All the variables with Kq magnitudes brighter than 
19.0 have been carefully checked on the reference (best 
seeing) image. The two Field 1 LPVs with periods around 
430 days, both of which are brighter than K s < 18.0 are 
located in a highly crowded areas. The fainter of the two 
is actually blended with a background galaxy and the 
brighter might be a blend of more than two stars. Also 
the third brightest Field 1 LPV could be too bright due 
to contributing light from the nearby stars. LPVs which 
arc blended with nearby stars should have smaller ampli- 
tudes due to larger relative contribution of the neighbour 
at minimum phases. Only one of the eight Field 2 bright 
LPVs (Kq < 19.0) is a possible blend of two stars. 



There are two sequences in the PL diagram in Fig. |SJ 
Multiple parallel sequences in the PL diagram of LPVs 
were reported by Wood & Sebo fiMfy . Wood QZMfy and 
Cioni et al. I)2001|l in the Magellanic Cloud. The more pop- 
ulated, longer period sequence corresponds to sequence 
C of Wood l|2000[l where Mira variables are found. The 
shorter period sequence B discussed by Wood ( 2000) is ap- 
proximately 1.3 mag brighter, while the shorter sequence 
in Fig. lies only 0.5-1 mag brighter. Moreover, it is is 
somewhat steeper, indicating a possible problem of alias- 
ing periods. 

Actually, a close inspection, shows that most of the 
stars on the shorter period sequence have periods close 
to 1/2 year and can be fitted similarly well with periods 
that are twice as long, and which place them on the Mira 
sequence. It is interesting to note that this increases also 
the amplitude of the fit. An example is shown in Fig. El In 
the upper panel the light curve is folded with the original, 
shorter, period and with the longer period on the bottom. 
Carefully phased additional observations of these variables 
would be necessary to determine their periods unambigu- 
ously. The period distribution of the LPVs and the com- 
parison between the two fields is discussed by Rejkuba et 
al. H2003bll . 

In Table all the variables that could be fitted with 2 
aliasing periods are listed. Columns 1 to 6 list: field identi- 
fication and ID number, original period, semi-amplitude, 
reduced x 2 °f the sine-curve fit, significance of the origi- 
nal period and the mean K s magnitude obtained from the 
sine-curve fit. The new period, semi-amplitude, reduced 
X 2 of the sine-curve fit with this new period and the as- 
sociated average K s magnitude, are listed in columns 7, 
8, 9 and 10, respectively. Additionally I have checked for 
possible aliasing periods all the LPVs that are located be- 
low the Mira PL relation. In most cases no shorter period 
could provide a satisfactory fit to the data, except for 2 
LPVs in Field 1 and 5 in Field 2. Fainter average iVband 
magnitudes for the LPVs that lie below the Mira PL rela- 
tion are due to extinction within their circumstellar shells 
(see also Kiss & Bedding 2003). In another 6 cases x 2 
of the sine-curve fit was improved and a new amplitude 
and mean magnitude were calculated, although the period 
remained very close to the original value. Apparently run- 
ning the sine-curve fitting algorithm in an automatic way 
failed to find the best fitting amplitudes in a few cases. 
These are also listed in Table |21 In the following distance 
determination these LPVs are not used. Fig. displays the 
PL diagram in which aliasing periods {P new ) from Tabled 
are used. 

A well-defined PL relation has been found for Miras in 
the Large Magellanic Cloud (Glass & Lloyd Evans Q981 
Wood UnSHJ), the Small Magellanic Cloud (Cioni et 
aL lSoTEty . the Galactic Bulge (Glass et al. lTMB) . the solar 
neighbourhood (van Leeuwen et al. I1997fl and in Galactic 
globular clusters (Feast et al. 2002 ) . The relation holds for 
both Mboi and Mr. Since Miras are very luminous, their 
tight PL relation makes them interesting for distance de- 
termination to other galaxies. 
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2.4 2.6 
log P (days) 

Fig. 8. Mira PL relation in NGC 5128 for variables with 
best determined periods and colors bluer than J—K < 1.4. 
The solid red line is our best fit to the Mira sequence. The 
dotted green line is the Feast et al. (2002) relation. 

Calibration of the if -band PL relation relies on the 
LMC PL relation for Miras. The Feast et al. ifE^ fit to 
the LMC Mira PL relation is: 



M K = -3.471ogP + /3 



(5) 



The zero point of (3 — 0.91 has been derived by Feast 
(1996) using LMC Miras. More recently Whitelock et 
al. I|2000|l used Hipparcos parallaxes for Solar neighbour- 
hood Miras to derive (3 = 0.84±0.14 and Feast et al. P3U2|> 
obtained (3 = 0.93 ± 0.14 for globular cluster Miras us- 
ing the cluster distances determined from subdwarf fit- 
ting with Hipparcos parallaxes of subdwarfs (Caretta et 
al. I2000[) . The straight mean of the two Hipparcos based 
zero points yields (3 = 0.88 ± 0.10 and a Large Magellanic 
Cloud distance modulus of 18.60 ± 0.10. Preferring in- 
stead a somewhat shorter LMC distance, (to — M)q AIC = 
18.50 ± 0.04 (Alves et al. IMI2) . I use a zero point of 
(3 — 0.98 ± 0.11 (where uncertainty in the LMC distance 
modulus has been added to the zero point error in quadra- 
ture). 

Glass & Lloyd Evans (2003) have re-reanalysed the 
Mira PL relation using the MACHO data for the same 
stars as Feast et al. (|1989ll in the LMC. They conclude that 
the periods of these variables did not change significantly 
and remained essentially constant over 2-3 decades. Their 
best fit to the if -band Mira PL relation is: 

if = -3.52(±0.21)logP + 19.64(±0.49) a = 0.13 (6) 

In order to determine the distance to NGC 5128 only 
Miras with colors J s — K s < 1.4 have been selected. This 
ensures that reddening from circumstellar shells that may 
be present in redder LPVs would not dim the if -band 
magnitudes and also matches the colors of Mira variables 
used to derive the Mira PL relation in the LMC ((J—K) < 
1.5, Feast et al. 1989). Moreover, only those LPVs that 



had the most regular periods, and could be fitted with a 
sine-curve with x 2 < 5 were chosen. They are plotted in 
Fig. |H1 together with the best linear least-square fit to the 
data (thick solid line): 

K = -3.37(±0.11)logP+28.67(±0.29) (N = 240) a = 0.20(7) 

The fits to Field 1 and Field 2 LPVs are identical within 
the errors. Largely, the error in the constant term is due 
to the extrapolation of the relation to the period of 1 day. 
The mean logP is ~ 2.5. If instead the fit is made at the 
mean period with the following functional form: 



K = A(\ogP-X) 



B 



(8) 



the error of the constant term drops to only 0.02. 

The slope of the PL relation in NGC 5128 agrees well 
with that in the LMC. It is within la of the slope of the 
Feast et al. (1989) relation and 1.5c of the Glass & Lloyd 
Evans l|2003|) PL relation. This is yet another piece of ev- 
idence that the slope of the PL relation in the if -band is 
universal, lending confidence to the distance determina- 
tion through the PL relation. 

In order to determine the distance modulus of 
NGC 5128 from the Mira PL relation, one has to assume 
that the slope is universal and use the same value as in 
the calibrating PL relation (i.e. that in the LMC). It is 
straightforward then to calculate the distance modulus 
from the zero point difference of the two PL relations. 
Restricting the slope to —3.47, and determining the zero- 
point of the relation at (logP — 2.3), approximately the 
mean period of the LMC stars used to determine the re- 
lation (Feast et al. lT58"9l Glass & Lloyd Evans the 
best fitting zero point is 28.94 ±0.03 (dotted green line in 
Fig.[SJ) with the same RMS of the fit. Using the above dis- 
cussed zero point of (3 = 0.98, I derive a distance modulus 
to NGC 5128 of (to - M) = 27.96 ± 0.11. If instead the 
slope of the PL relation is fixed to —3.52 (Glass & Lloyd 
Evans l27)03[) . the resulting distance modulus to NGC 5128 
is (to-M)o = 27.93 ±0.11. 

mm 



Feast et al. jl989|l as well as Glass & Lloyd Evans 
(2003) noted the group of stars with periods in excess of 
~ 420 days which lie some ^0.7 mag above the PL rela- 
tion. Due to the fact that a few of them have been found 
to be Li-rich, it was suggested that they are more mas- 
sive, hence younger and are in the Hot-Bottom Burning 
phase (Smith et al. H995fl . So, restricting the LPV data in 
NGC 5128 to LPVs with periods shorter than 400 days, 
where the above mentioned LMC PL relations have been 
calibrated, distance moduli of 27.95±0.12 and 27.92±0.12 
are derived adopting slopes of —3.47 and —3.52, respec- 
tively. These distance moduli are virtually indistinguish- 
able from that determined from the full data set and they 
point out that most of the long period stars do follow the 
same Mira PL relation as shorter period Miras. They are 
not overluminous and in the hot-bottom burning phase. 

There are a few dozen stars in Fig. [5] that are J; 2a 
above the PL relation. Some of them might be brighter due 
to contributing light from their neighbours, but some show 
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Table 3. Summary of the distance determinations in this 
work. 



DM (mag) 


Method 


Field 


27.87 ±0.16 


A"-band RGB tip 


Fl 


27.9 ± 0.2 


//-band RGB tip 


Fl 


27.96 ± 0.11 


/('-band Mira PL relation 


Fl & F2 



indications of humps in their light-curves. Similar humps 
have been found in MACHO light curves of the LMC Li- 
rich stars by Glass & Lloyd Evans (2003). However, before 
being able to place firm conclusions, accurate light curves 
over several periods should be determined for a larger sta- 
tistical sample of Li-rich hot-bottom burning AGB stars. 

5. Conclusions 

The tip of the RGB was detected in Field 1 at K s = 
21.24 ± 0.05, yielding a distance modulus of NGC 5128 of 
(m - M) = 27.87 ± 0.16. The comparison of the //-band 
RGB tip luminosity in the Galactic Bulge and NGC 5128 
implies a similar distance modulus (Table |3J). The RGB 
tip in Field 2 is not a sharp feature due to a brighter com- 
pleteness limit, larger photometric errors and the presence 
of blends and AGB stars. 

A large population of stars above the tip of the RGB 
contains 2176 stars in the outer halo field (Field 1) and 
6072 stars in the inner halo field (Field 2). Subtracting 
these foreground sources, detected LPVs, as well as max- 
imum probable number of blends of two RGB tip stars, 
there are some 1150 and 150 non- variable stars brighter 
than the first ascent giant branch tip in the two fields. 
LPVs account for 26% and 70% of the AGB popula- 
tion in Fields 1 and 2, respectively. The high luminos- 
ity [M-k < —8.7) achieved by AGB stars is a sign of an 
intermediate-age population. 

The first Mira period-luminosity relation outside the 
Local Group is presented. Miras with the best determined 
and most regular periods that do not have red colors 
(J s — K s < 1.4) were used to determine the distance of 
NGC 5128 from a period-luminosity relation. I derive the 
distance modulus of 27.96 ±0.11, adopting the LMC dis- 
tance modulus of 18.50 (Alves et al. 12002(1 . 

The mean distance of the Miras, obtained from fields 
on both sides of the center of NGC 5128, is slightly larger 
than the mean distance of the red giants in Field 1. This 
indicates that the orientation of NGC 5128 is such that the 
north-eastern part of the galaxy is closer to us. However, 
the relative distance of the two fields is smaller than ~ 
0.03 Mpc and negligible compared to the distance from 
the Milky Way. The mean of the two methods yields a 
distance to NGC 5128 of (m - M) = 27.92 ± 0.19 {D = 
3.84±0.35 Mpc), very close to the (m -M) Q = 27.98±0.14 
measurement of Harris et al. I|1999H and within la of the 
distance modulus 27.8 ± 0.2 result of Soria et al. (1996). 
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Table 2. LPVs with possible aliasing periods. Original field identification and ID number, period, semi-amplitude, 
reduced x 2 of the sine-curve fit, significance of the original period and the mean K s magnitude are listed in columns 
1 to 6. The new, aliasing, period, semi-amplitude, reduced \ 2 of the sine-curve fit and the associated average K s 
magnitude, are listed in columns 7, 8, 9 and 10, respectively. 
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